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1. Introduction 

Various graph coloring problems have been considered in the last two decades in 
relation with statistical mechanics. On regular lattices, many of these have proved to 
be integrable and exactly solvable ffl. Among these is the classical tri-coloring problem 
of the bonds of the triangular lattice, solved by Baxter in [Q: in how many ways can 
one color the lattice bonds with three colors so that any triangle has its three edges of 
a different color. Recently a remarkable connection was made between this problem and 
that of folding of the triangular lattice ||, counting the different ways to fold the lattice 
onto itself along its bonds, a toy model for the crumpling of two-dimensional polymerized 
membranes. A natural extension of these problems consists in their generalization to 
random triangulations. This may be viewed as a coupling of the model to two-dimensional 
quantum gravity, or also as a way to describe fluid membranes. The bond tri-coloring of 
random triangulations has first been addressed in || and solved in |5| for a restricted model 
where the colors always occur in the same cyclic order for all triangles. The complete bond 
tri-coloring problem without restriction was solved very recently in || . The connection to 
the folding of random triangulations requires unfortunately an additional constraint: the 
number of triangles around each vertex must be even. The folding of random triangulations 
remains an open problem. 

In this paper, we address a different kind of (bi- or tri-)coloring problems, where we 
now color the vertices or the faces of some random graphs, representing tessellated surfaces 
of arbitrary genus. Starting from the vertex bi-coloring of arbitrary graphs, our main result 
will eventually be rephrased as solving the tri-coloring problem of the vertices of random 
triangulations of surfaces. 

Our starting point is the computation of the following integral over two Hermitian 
matrices M and R of size N x N (see for a review) 

Z{p,q,g;N) = J dMdRe~ NTlV ( M ^ p ' q ' 9 ) (1.1) 
with the non-polynomial potential 

V(M,R;p,q,g) = pLog(l — M) + <?Log(l — R) + gMR (1.2) 
and with respect to the Haar measure over Hermitian matrices H 

N 

dH = c N (g)Y[dHii ] \ dReH tJ r fflmify (1.3) 

i=l l<i<j<N 
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normalized in such a way that Z(0, 0, g; N) = 1 for all g and N. The parameters p, q, g are 
a priori any real numbers. This model, with its logarithmic potential, can be viewed as a 
two-matrix generalization of the one-matrix Penner model H. For p = q, it has already 
been studied in |J in connection with the Kazakov-Migdal model ||10|| . 

Upon setting g = 1/t and rescaling M — > yiM, R — > y/tR in the integral (|1.1|) , we 
may rewrite ( |1 . 1|) in the same form, but with the potential 

W(M,R;p,q,l/t) = pLog(l - VtM) + gLog(l - VtR) +MR (1.4) 

and with a different normalization of the measure, still ensuring that Z(0, 0, 1/t; N) = 1. 
Note that this last condition is equivalent to demanding that lim^o Z(p, q, 1/t; N) = 1. 
If we think of t as a small parameter, we may at least formally expand the integral Z in 
powers of t. This expansion is expressible as a sum over Feynmann (fat)graphs, made of 
double-lines oriented in opposite directions and carrying a matrix index i e {1,2, ...,N}, 
involving two types of vertices (M-vertices weighed by Np and /^-vertices weighed by Nq, 
irrespectively of the valency) connected by one type of propagator ((MR) propagator, 
weighed by t/N). As usual in matrix integrals, each graph also receives a contribution N 
per oriented loop, from the summation over the running matrix indices. Hence we finally 
get the expansion 

^ t^(r) N x(r) v v p v q 

z^Vt,N)= E | Aut( r)| ^ 

fatgraphs T 1 K n 

where E(Y) denotes the number of edges of T, x(r) its Euler characteristic, V p , V q its 
numbers of M and /^-vertices respectively, and |Aut(T)| the order of its symmetry group. 
The sum may be restricted to connected fatgraphs by considering Log Z instead. 

If we choose to color the vertices of the fatgraphs with two different colors, say white 
for the M-vertices and black for the -R-vertices, the occurrence of only mixed propagators 
(MR) implies that the graph is bi-colored, namely two adjacent vertices have different 
colors. The white (resp. black) faces receive the weight p (resp. q). We can also consider 
( |1.5| ) as an expansion over the dual fatgraphs, with now black and white alternating faces. 
So we may finally view Log Z as the partition function for connected bi-colored fatgraphs 
of arbitrary genus. 

When compared to the class of exactly computable multi-matrix integrals, ( |1 . 1| ) is 
readily seen to be reducible to an integral over the eigenvalues of the matrices M and R, by 
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use of the Itzykson-Zuber integral formula [[11]. However, we expect the non-polynomiality 
of the potential (|1 . 2| ) to translate into some interesting non-locality of the model, which 
makes its solution slightly more subtle than that of the polynomial case. 

The paper is organized as follows. The model ( |ITT| ) is solved in Sect. 2, using standard 
orthogonal polynomial techniques. The solution involves the computation of partition 
functions Z n with the same definition as ( |1 . 1|) , but where the integration extends over 
n x n Hermitian matrices, for all values of n between 1 and N. The varying size n of 
the matrices introduces a new scaling variable z = n/N in the problem. Remarkably, 
the solution exhibits an explicit (p, q, z) symmetry relating this size scaling variable z to 
the coupling constants p and q of the model. This symmetry is explained by identifying 
our problem with that of tri-coloring the vertices of random triangulations, with the third 
color (say grey) receiving precisely the weight factor z. The solution is further explored 
by taking the large iV limit, in which only planar (genus zero) diagrams contribute. The 
result is a non-linear partial differential equation for the free energy, involving derivatives 
wrt t and z. To solve it, we first derive initial data when z —>■ 0, in which case the free 
energy may be explicitly computed (the detailed derivation is done in Appendix A). In 
Sect. 3, we present an alternative solution of our model, by first writing the integral fll.lj ) 
as a determinant, and then deriving a (discrete Hirota) bilinear equation which determines 
the partition function completely. In the large iV limit, this takes the form of yet another 
non-linear partial differential equation for the free energy, involving derivatives wrt p, q, 
z and t. Combining the results of Sects. 2 and 3., we are able to compute explicitly the 
free energy in the case p = q = z. In Sect. 4, we follow a third route to express the 
partition function Z n of our model as an explicit multiple sum over strictly increasing 
sequences of integers. We further extract the large iV behavior of these sums by applying 
the saddle-point approximation, which yields yet another expression for the planar free 
energy. The latter permits to investigate the critical properties of the model. We find a 
first order critical point whenever p, q and z have the same sign, with a generic singularity 
Log Z ~ (t*— t) 5 / 2 , characteristic of the universality class of pure two-dimensional quantum 
gravity. The all-genus double-scaling limit of the free energy is also derived, and shown 
to be governed by the Painleve I differential equation. We conclude in Sect. 5 with a 
discussion on the connection of our model to folding problems and with some remarks on 
the emergence of discrete Hirota bilinear equations in multi-matrix models. 
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2. Orthogonal Polynomial Solution 

2.1. Preliminaries 

As mentioned above, we first reduce the integral ( |1 . 1|) to an integral over the eigen- 
values of M and R, denoted by rrii and respectively. We simply perform the change 
of variables (M;R) — ► (m,U;r,V), with M = UmU^ and R = VrV\ U,V two unitary 
matrices, and m =diag(rai, mjv), r =diag(rx, rjy)- The Jacobian of the change of 
variables is simply A(m) 2 A(r) 2 , where A(m) = rii<j( m * ~ m j) denotes the Vandermonde 
determinant of the matrix m. The integral over the unitary matrices U and V only involves 
the cross-term —NgTt(MR) in the potential and reads 

r det \e~ N9mir A 

/ dUdVe~ N ^ MR ) = d N (g) 1 , A T J ~ N ( 2J ) 
J A(m)A(r) 

by direct use of the Itzykson-Zuber integral formula. The precise value of the normalization 
constant d^(g) is not important here. Substituting this into the integral (|1.1|), and using 
the antisymmetry of the Vandermonde determinant, we finally get the reduced integral 

Z(p,q,g;N) = J dmdrA(m) A^e"™^™™'^ (2.2) 

where the integral extends over the 2 N real variables m, , r j , and the iV-dimensional mea- 
sure over diagonal matrices I =diag(/i, ...Jn) 

N 

dl = Ve N (g) 11 dk (2.3) 

is normalized so as to ensure that Z(0, 0, g] N) = 1. The normalization factor ejv(fl') will 
be computed later. 

As in the standard orthogonal polynomial technique, we introduce two sets of monic 
polynomials p n (x) = x n -|-lower degree, and p n (y) = y n +lower degree, for n = 0, 1, 2, 
which are orthogonal with respect to the one-dimensional measure inherited from ( |2.2| ), 
namely 

{Pn,Pm) = J dxdye- NV ^' y ' p ^p n (x)p m (y) = h n S m , n (2.4) 

for some normalization factors h n (p, q, g; N), entirely fixed by the orthogonality condition. 
Using the multi-linearity of the determinants, we may rewrite 



A(m) = det[m^ 1 ]i<ij<jv = det [pj-i(mi)]i<i,j<jv 

(2.5) 



A(r) = det[r^ ]i<ij<jv = det \pj-i(ri)]i<i,j< 



N 
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and expand the latter. Using the orthogonality relations, the partition function is finally 
rewritten as 



jV-l 



Z(p,q,g;N) = N\e N (g) ] [ K 



i=0 



(2.6) 



2.2. P, Q operators and relations 



Next we write recursion relations linking the h^s by introducing the operators Q and 
Q of multiplication by an eigenvalue of respectively M and R, namely 



(QPn)(x) = XPn(x) (QPn)(y) = VPniv) 



(2.7) 



and the operators P and P defined by 

(PPn)(x) 



X 



dp n {y) 



(2.! 



Note that these differ slightly from the usual definitions for the polynomial potential case 
(where P and P are simply the derivatives wrt the corresponding eigenvalue), and that 
they satisfy different (non-canonical) commutation relations 



[P,Q] = Q-l [P,Q] = Q-l 



(2.9) 



Equations for the h^s are derived by considering various matrix elements of the op- 
erators P and P. 



(Pp ni p m ) = -(p n ,Pm) + N((x - l)d x V(x,y;p, q 1 g)p ni p m ) 
= (Np-l)(p n ,p m ) + Ng((Q-l)p n ,Qp m ) 



(2.10) 



where we have performed an integration by parts wrt x. Note that the particular definition 
of P is ad-hoc to eliminate all non-polynomial dependence on Q in the final equation. 
Denoting by the adjoint of an operator A wrt the bilinear form (/, Ag) = (A^ /, g), we 
can recast (|2.10[) into an operator identity 



Analogously, we have 



P+(1-Np) = NgQ\Q - 1) 



P+(1-Nq) = NgQ\Q-l) 



(2.11) 



2-12) 
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The operator relations (|2.9|), ( |2.11| ) and ( |2.12|) determine the h^s entirely. 

Note that in our definitions, P, Q and P, Q act on two different bases (p n ) and (p n ) 
of the space of polynomials of one variable. Let us now write the actions of the operators 
P, Q, and on the same basis (p n ). We have 



(QPn)(x) 


^ Qn,kPn—k(x 




k>-l 


(PPn)(x) 


= ^2P n ,kPn-k(x) 




k>0 


{Q ] Pn){x) 


= ^2 Qn,kPn+k(x 




k>-l 


{P ] Pn){x) 


= S ^Pn,kPn+k{x) 




k>0 



where the first terms are given by 



Qn,-1 
Pn,0 

Qn.-l 

Pn,0 



n 



h, 



h n -i 
n 



(2.13) 



;2.i4) 



The first two lines of (|2.14|) follow directly from the definitions (|2.8|) , whereas the last two 
lines are easily derived from the identities 



(&Pn,Pn-l) = (Pn,QPn-l) = h n = Qn-xhn-x 

(P^p n ,p n ) = (p n ,Pp n ) = nh n = P n ,oh n 



(2.15) 



We obtain all the relations between the coefficients in ( |2.13| ), by expressing ( |2.9| ) and 
(pH|)(|2TT2|) as 



[P,Q] 


= Q-l 




[Q^P j ] 


= Q f -l 




P + (1 - Np) 


= NgQ\Q 


-1) 


P f + (1 - Nq) 


= Ng(&- 


i)Q 



;2.i6) 



and letting these operators act on (p n ), namely 

Qn,k &k,0 ^ ^ Qn,mPn—m,k — m Pn,k — mQn — k-\-m 

m>-l 

Qn,k $k,0 ^ Pn,k—mQn-\-k — m,m Qn,mPn-\-m,k- 



(2.17) 



m> — 1 



and 



;i - Np)6 k , + P n , k + NgQ n> - k = Ng ^ Q 

m> — 1 

Pn, — k 

+ NgQ n , k = Ng Q 

n — m,m — kQn,r 



(2.18) 



m> — 1 

This gives an infinite set of equations for the coefficients of (|2.13|) . A drastic simplification 



of ( |2.18| ) follows from subtracting the third from the fourth line of ( |2.1(j| ), cancelling the 
crossed term Q^Q. This results in the operator identity 

A = P-NgQ-Np = (P - NgQ - Nq) j = I f (2.19) 

where we have formed two convenient linear combinations A and A of the P, Q and P, Q 
respectively. When expressed on the basis (p n ), ( |2.19|) yields 

Pn,k-NgQ n , k -Np5 kfi = P ni _ fe - NgQ n> _ k - Nq5 kfi (2.20) 

Taking k > 1 (resp. k < —1) in this formula, in which case the rhs (resp. the lhs) vanishes, 
this gives 

P n ,k = NgQ njk P n , k = NgQ Uyk (2.21) 
for all k > 1. Taking now k = 1 and k = — 1 in ( |2.20| ) leads to 

Pn,x-NgQ n ,x = -NgQ n -! 

(2.22) 

Pn.x — Ng Q n x = -NgQ ni -t 



while k = 0, ( ggjj) gives 



Qn,0 — Qn,0 — ~ (2.23) 

9 



But even after eliminating P n , k , Q n ,k and P n ,k, ( |2.17| )( pTT8| ) remain an infinite set of 
quadratic equations for the Q n ,k- 



In the next section, we shall derive a set of equations involving only the first values of 
the index k, at the cost of taking derivatives wrt g. We shall make use only of Eqs. ( |2.14| ) 
and ( |2.23| ) , together with the first line of ( |2.17| ) for k = 0, which yields 

Qn,0 — 1 = Pn+1,1 — Pn,l (2.24) 

This equation enables us to solve for P nj i in terms of the Qi,o, given that Po,i = by 
definition. Introducing the quantity 

n-1 

J2Qk,o (2.25) 



k=0 



we find that 



2.3. The main recursion relation 



Pn,i = v n -n (2.26) 



In this section, we take advantage of the simple dependence of Z on g to derive 
differential equations for the first coefficients Q n ,k, which, supplemented to the results of 
the previous section, yield a recursion relation for Z. 

Let us evaluate the following derivatives wrt g (the notation d g stands for d/dg) 

dg(pniPn-l) = = (dgPn.pn-!) + (pn^dgpn-!) - N(Qp n ,Qp n - l ) 

= {d g p n ,p n - X )-N{QtQp n ,p n -. x ) (2.27) 
d g {p n ,p n ) = d g h n = -N(Q ] Qp n ,p n ) 

where we have used the fact that d g p n has degree < n — 1, and the orthogonality of the 
p's and p's. To compute (d g p n ,p n -i), let us write p n (x) = x n — A n a; n_1 + ... Using the 
definition of Q ( |2.7|) , we easily see that A n = v n , given by (|2.25|) . Indeed, we write 



Pn+l(x) = X n+1 -X n+1 X n +0(x n - 1 ) 

= QPn(x) - Q nfi x n + Oix"- 1 ) (2.28) 
= x n+1 -(A n + Q n ,o)x n + 0(x ri - 1 ) 

and the result follows. Hence 

{dgp n ,p n -i) = -h n _ x dgV n (2.29) 



8 



Finally, let us express N(Q^Qp n: p m ) by using ( |2.11|) , namely that NQ^Q = NQ^ + (P + 
1 — Np)/g, hence 



N(&QPn,Pn] 



/i n _i(A^g n ,_i + -p n>1 ) 

K(NQ n , + -(Pn 1 o + l-Np)) 



(2.30) 



and finally, using the explicit values ( |2.14| ) and ( |2.26| ), we get the differential equations 



(d g + -)v n 
9 

d g Log h n 



9 h n -x 
1 



-(Np-n-l)-NQ n , 

= -(Nq-n-l)-NQ ny0 
9 

where we have used ( |2.23| ) to get the last line. Let us introduce 

n-l 

F n = ^Log/i fc 

k=0 

so that 

Log Z(p,q,g;N) = F N + Log(e N (g)N\) 
then after summation, the second equation of ( |2.31| ) becomes 



8 F 



1 

2^ 



n(2Nq -n-l) - Nv r 



After elimination of v n between (|2.31| ) and (|2.33|) , we are left with 

n 

Ng 



-^(d g + ~)d 9 F n 



h n n 



h 



n-l 



In terms of the variable t = 1/g, and upon the definition 

n = 1,2,... 



h n -i 
a = h 

we finally have the system of equations 



n-l 






- k)Loga k 


k=0 




t 

Oi n = n— + 


^ t fF n 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



^2.37) 
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Given the initial data ao = ho(p, q,l/t; N), this set of equation gives a recursion relation 
which determines F n for all n, and in particular F^. 

The normalization factor e^ig) in ( |2.33| ) still has to be fixed. It is obtained from the 
solution (an \ Fn°^) of the system (|2.37|) for the case p = q = 0. In that case, the initial 



data reads 



We then have 



for n > 1, and 



4°) = h o (0,0,l/t;N) = J e- Nxy ' t dxdy = (2.38) 



a<P = n-L (2.39) 



4 0) = L °s4 + £(^) L °§(S) 

( / + \ n(n+l)/2 n-1 



(2.40) 



3=0 

Now the normalization e^ig) in ( 2.33 ) ensures that Z(0, 0, AT) = 1, hence we have 



Log Z(p,q,l/t;N) = F N - F$ ] (2.41) 
Upon comparison with ( |2.33|) , we deduce that 



1 

(t/N) N(N ^ /2 UU^ 



e "(V*) = „,_^ rilWo -v - (2.42) 



Introducing the renormalized quantities 



a n = for n = 1,2, ... 

nt 

Na 



(2.43) 



and the free energy 



such that 



/„ = F n - (2.44) 



f N = Log Z(p, q,l/t;N) (2.45) 
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we finally have the recursive system, for n > 1 




(2.46) 



with the initial data 



a 



N 
T 

£ 

fc>0 



y cfedy(l - ^-^(1 - y )- N i e - Nxy/t 

i r(JVp + fc)r(iVg + k) f t 

kl T(Np)T(Nq) 



(2.47) 



Note that if — Np and/or — Nq is a positive integer, the above sum truncates. 

At this stage, it is interesting to note that the quantity f n which appears in the 
recursion relation ( |2.46|) has the following simple interpretation. Indeed, for n > 1, we 
have 

/„ = LogZ n (2.48) 



where Z n = n\e n (g) YYi=o is the integral (TTT) taken over n x n matrices, but with the 
same potential (|1.2| ), with iV as a prefactor. The quantity f n is thus the corresponding 
free energy. The normalization factor e n (g), which ensures that Z n = 1 for p = q = is a 



straightforward extension of (|2.42 ) 



e n (g) 



(Ng) 



n(n+l)/2 



(2.49) 



In terms of fatgraphs, we now have a weight n instead of N per oriented loop, leading to 

(n)^) ^(F) N x(T) p V pq V q 



Z n (p,q,l/t;N) = J2 

fatgraphs V 



|Aut(T)| 



(2.50) 



where F(T) denotes the number of faces of the graph. Therefore we may think of the new 
scaling variable z = n/N as an activity per face of the bi-colored graph. 
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2.4- Solution of the main recursion relation 



To solve (|2.46| ), we first eliminate a n by noting that 

Z n +lZ n -\ 



a, 



Zl 



(2.51) 



for all n > by setting Z_i = Zq = 1. This is obtained by inverting the relations 
=0 a fc • I n terms of Z n , the recursion relation (|2.4d| ) reads 



Zn-lZn+l 



nN 



{td t ) 2 hogZ r 



(2.52) 



for all n > 1. The initial data are Z_i = Zq = 1, and Zi 
Let us first solve (|2.46|) in the case 



ao as in (|2.47|) 



a = Np = -1 



(2.53) 



and b = Nq arbitrary. The initial data ( |2.47 ) reads 



ao = 1 — 6 



N 



(2.54) 



The solution of ( |2.52j ) takes the form 

Z " = l + E(^)6(6-l)...(6-* + l)(~)' 



(2.55) 



for n > and Z_i = 1, as easily checked by explicitly differentiating the rhslil. In particular, 
the solution for n = N reads 

Z(-l/N,q,l/t-N) = l + J2(jyNq)(Nq-l)...(Nq-k + l)(-^J (2.56) 



A similar reasoning permits to solve the case 



a = Np = 1 



(2.57) 



1 This result will also be proved directly in Sect. 5. 2 below. 
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and b = Nq arbitrary, with now 



i + £(' B + J" 1 W + i)--(»+*-i-" 



k=l 



N 



(2.58) 



also satisfying ( |2.52| ). The corresponding solution for n = N reads 



Z(l/N,q,l/t;N) = l + f:{ N + k k ~ 1 ) 



N + * 1 ^Nq(Nq+l)...(Nq + k-l)^) (2.59) 



More generally, we can use the recursion ( |2.52[ ) to expand f n = Log Z n (p, q, l/t\N) 
as a formal power series of t 



f n (p,q,l/t;N) = ^w n , fc (a,o) (-^j 
k>i ^ ' 



(2.60) 



with a = Np, b = Nq. The initial conditions Zq = 1 and Z\ = ao translate into the 
conditions 



uj ,k(a, b) = 



fc>i 



(2.61) 



with ao as in ( f2.47] ). The recursion relation ( ^.52| ) turns into a recursion relation for the 



<^n+l,fc — 2w nj fc — U>n-l,k 
(~ I) 7 "" 1 



r>l 



fc 1 , .... few. > 1 



(2.62) 



Together with the initial conditions ( [2.61| ), this determines the coefficients uj n ^{a,^ b) com- 
pletely. It is easy to show that these are polynomials of a, b, and n for each value of k. 
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The first few of them read 



a; n ,i = nab 
nab 



W n ,5 



[n + a + b) 
nab 

= —(n 2 +3(a + b)n + a 2 + 3ab + b 2 + l) 

O 
71 d u 

= — (n 3 + 6(a + b)n 2 + (6a 2 + llab + 6b 2 + 5)n 

+ (a + b)(a 2 + 5ab + b 2 + 5)) 
nab 



n 4 + 10(a + o)n 3 + 5(4a 2 + lla6 + 46 2 + 3)n z 
u 

+ 5(a + b) (2a 2 + 9ab + 2b 2 + 8)n 

+ a 4 + 10a 3 6 + 20a 2 6 2 + 10a6 3 + b 4 + 15a 2 + 40a6 + 156 2 + 



(2.63) 




(a) (b) 

Fig. 1: The construction of the tri-colored graph V (Jo) corresponding to the 
bi-colored graph T (a). Each face is subdivided into triangles, by adding a 
new central vertex, colored in grey. The face-weight z is re-affected to each 
such new grey vertex. 

Remarkably, the result ( |2.63|) reveals that each coefficient u> n ^{cL, b) is totally symmet- 
ric in the three variables (a, 6, n), or equivalently the three variables (p, g, z = n/N). Let us 
prove this property directly on f n (p, q,l/t; N) ( |2.60|) by considering its fatgraph expansion 
( 2.5 0|) (restricted to connected graphs). For each fatgraph F with black and white vertices 
in the expansion ( |2.50| ), let us construct one grey vertex in the middle of each face (see 
Fig.l), and connect it to each (black or white) vertex around the face. As black and white 
vertices alternate around each face, this creates an even number of edges. Each grey vertex 
receives a contribution z, whereas white vertices are weighed by p and black ones by q. 
We end up with a fatgraph V whose vertices are tri-colored (in white, black and grey), and 
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whose faces are all triangles with a vertex of each color. The automorphism group remains 
unchanged (|Aut(T)| = |Aut(T)|), so as the Euler characteristic (x(T) = x(r)), hence we 
get the expansion 

F Q V b z Vg N X 

f n (p, q ,l/t;N) = * t P lAt(f]l ( 2 - 64 ) 

conn, tricol. I \ / I 

fatgraphs T 

where F denotes the total number of (triangular) faces of T (F = 2E(T)), and V w , Vb, V g 
respectively its numbers of white, black and grey vertices. The form ( |2.64j) is explicitly 
symmetric in the three variables (p, q, z\ that is in the three variables (a, b, n). 

The a <-> b symmetry was clear from the integral ( |1 . 1[ ) . Another direct proof of the 
symmetry n <-* a will be given in Sect. 4. 2 below. As an additional check, we may compare 
Zi(a = n,b) = ao(a = n, 6, t) as given by (|2.47|) , to Z n (l, b) found in (|2.58|) : the two sums 
are indeed identical. 

2.5. Large N limit 

In the large iV limit, the free energy Log Z is dominated by the planar (genus zero) 
fatgraphs, hence a leading contribution 

Log Z(p,q,l/t;N) ~ N 2 F (p,q,t) (2.65) 

when — > oo. Introducing the rescaled variable z = n/N, we may assume that a n tend to 
a smooth function a(z;p, q, t) of z G [0, 1] in the large N limit. This is in agreement with 
the assumption that f n /N 2 also tends to a smooth function f(z;p,q,t). Indeed, ( |2.46j ) 
then implies that 

a(z;p,q,t) = 1 + ^(td t ) 2 f(z;p, q, t) (2.66) 

According to (|2.50| ), the continuous quantity z is nothing but a fugacity associated with 
the faces of the (planar) fatgraphs. 

The recursion relation ( |2.52| ) reads, for e = 1/N and f(z) = f(z;p, q,t) 

e N*(f( Z+ e)+f(z-e)-2f(z))+0(l) _ 1 = ±_ ^ 2 (JV 2 /(z) + Oil) ) (2.67) 

nN 

hence in the limit N ^ oo 



d 2 J = Log(l + ^(^) 2 /) 
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(2.68) 



The partial differential equation (|2.68| ) however does not determine f(z;p, q, t) completely. 
We need to supplement it by some initial condition. This is easily done by use of the 
formal small t series expansion of /. From the large N limit of ( |2.60|)([2.63|) (with a = Np, 
b = Nq) , we indeed find that 

f n = n ab^- + 0(t 2 ) 

N 1 ' (2.69) 

= N 2 pqzt + 0(t 2 ) 

hence we get the initial condition 

f(z;p,q,t) = pqzt + 0(t 2 ) (2.70) 
This can be used to feed the differential equation ( |2.68| ). Indeed, writing 

f(z;p,q,t) = Y^^z-p^q)^ (2.71) 

fe>l 



with ui = pqz, the equation (|2.68|) turns into a recursion relation 

d ^ = E kl..k 2 r u kl ...u; kr (2.72) 

7"^l ^ 1 J ■ ■ ■ : k r > 1 

T.k i = k-r 

for k > 2. This is nothing but the continuum limit of the recursion relation ( |2.62|) . 
Eq.( [2.72l) however involves two integrations wrt z at each step, hence the introduction of 
many integration constants. To fix them, we need to know the values of f(0;p,q,t) and 
d z f(0;p, q, t). The latter are readily obtained by analyzing the large iV limit of the initial 
condition ( |2.47| ), namely we write 



= x - 1 T(Np + k)T(Nq + k) ( t \ k 
1 ~ G ° ~ T(Np)T(Nq) \Nj 



k>0 

N I dae NS{a ' p ' g ' t) 



(2.73) 



o 



where we have replaced the summation over k = aN by an integral over the positive real 
numbers a, and used the Stirling formula to derive the effective action 

S(a;p,q,t) = (p + a)(Log(p + a) - 1) + (q + a)(Log(<? + a) - 1) + aLogt 

(2.74) 

— a(Log a — 1) — p(Logp — 1) — g(Log q — 1) 
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The integral is dominated by the saddle-point dS/da(a*) = 0, namely 



t(a* + p)(a* + q) = a* (2.75) 
We only retain the solution leading to a maximum of S, namely 

a * = 1 - \ + q) - h ^ ~ t{p + q))2 ~ 4pqt2 (2 - 76) 

Finally, the free energy fi = Log Z\ behaves for large N as 
fi(p,q,t) ~ AS(a*) 

= A(pLog(l + — )+ 3 Log(l + —)-«*) (2.77) 
= N(p t (p, q,t) 

Let us now interpret the result fl2.77f) . As n = 1, we have z = 1/iV. The fact that /j is of 
order N = N 2 z shows that 

f(z;p,q,t) = zvi{p,q,t) + 0{z 2 ) (2.78) 

for small z, hence that 

f(0;p,q,t) = 

d z f(0;p,q,t) = ipxip^t) = g)f* (2 ' 79) 

fc>i 

with (pi as in ( |2.76| )( pT77|) . The coefficients Hk(p, q) are computed in Appendix A, and read 

= F£(i)G-i) pV+w (2 - 8o) 

The quantity (jp, q, t) can be viewed as the free energy associated with bi-colored graphs 
with exactly one vertex, i.e. bi-colored trees. As clear from (|2.76|) , it develops a singularity 
at a critical value for t equal to 

U(p,q,0) = ~ -2 (2-81) 

For t — > t* (p, q, 0) , we find (see Appendix A) 

^i(p,g,t)~(t*-*) 3/2 (2.82) 
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As we shall see however, this behavior for the free energy, valid only in the limit z — > 0, is 
not generic. 

With the initial data ( |2.79| ) ( pT80|) , we may now write the differential equation ( |2.72j ) 
in integrated form with u>k(z) = Uk(z;p, q) 

rz i>x / l) r— 1 

J ° J ° r>l *i,...,kr>l 

£fc i= fc-r 

This is now a true recursion relation on k, which determines all u/s. With the initial 
condition u Q (z;p, q) = 0, we find 

u>\ = pqz 

pqz 

uj 2 = -Y{z+p + q) 

^3 = ^-|-(z 2 + 3(p + q)z + p 2 + 3pq + q 2 ) 

U4 = ?-^(z 3 + 6(p + q)z 2 + (6p 2 + llpq + 6q 2 )z + (p + q){p 2 + 5pq + q 2 )) 

U5 = ^( z 4 + io(p + g ) 2 3 + (20p 2 + 55^ + 20Q 2 )2 2 
5 

+ 5(p + q)(2p 2 + 9pq + 2q 2 )z + p 4 + 10p 3 q + 20p 2 q 2 + 10pq 3 + q 4 ) 

U6 = ^(z 5 + 15(p + q)z 4 + 5(10p 2 + 27pq + 10q 2 )z 3 
6 

+ 2(p + q)(25p 2 + 106pq + 25q 2 )z 2 + (15p 4 + 135p 3 q + 262p 2 q 2 + 135pq 3 + 15q 4 )z 
+ (p + q)(p 4 + 14p 3 q + 36p 2 q 2 + Upq 3 + q 4 )) 

(2.84) 

For z = 1, this gives the first few terms in the expansion of the planar free energy ( |2.65| ) 
as a formal power series of t. In this expansion, the coefficient of t k p r q s corresponds to the 
planar connected (dual) black and white fatgraphs with r white faces, s black faces, and k 
edges separating them, weighed by the inverse of the order of their symmetry group. For 
z arbitrary, the coefficient of z m singles out those black and white graphs with exactly m 
vertices. Since we are dealing with planar graphs, one has m + r + s = k + 2. This explains 
why Uk(z;p,q) is a homogeneous polynomial of p, q, z, of degree k + 2. Note also that 
(|2.84|) agrees with the large iV limit of (|2.63| ), i.e. 

Jim -r^^zN,k(Np,Nq) = u k (z;p,q) (2.85) 

iv — >oo iv 

Finally, the result Q2.84 ) is explicitly symmetric in the three variables (p, q, z), in agreement 
with the expected symmetry 

f{z;p,q,t) = f(z;q,p,t) = f{p;z,q,t) (2.86) 
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3. Determinant form and discrete Hirota equation 

3.1. The partition function as a determinant 

Throughout this section, we use the shorthand notation Z n (a,b) for the partition 
function Z n (p, q, g; N) (|1 . 1|) for an integral over n x n matrices, and with a — Np, b = Nq. 

We start from the reduced integral (|2.2|) (|2.3|), extending over the 2n eigenvalues m; 
and Ti 

/n 
J] drmdnil - m,)- a (l - ri )- b A(m)A(r)e- N9 ^ min (3.1) 
i=i 

From the definition of the Vandermonde determinant, we have A(l-m) = (-l) n ( n - 1 )/ 2 A(m), 
hence we may replace m — > 1—m and r — > 1— r in both Vandermonde determinants without 
altering Q3.1|) . Moreover, 



\[{l-m i )- a A{l-m) 



det 



i=i 



(1 - Wi) 



j-a-1 



- l=*)i= n 

W(i)- a -l 



(3.2) 



and we have an analogous formula for r (with b instead of a) . Substituting these expansions 
into (PTT|) , we get 

Z n (a,b) = e n (g) ^ sgn(£7r) 



<T,T£S r , 



- N gmiVi 



x \\ / dTTiidriCl-^^-o-i^.^rW-ft-ie ......... (3 3) 

1=1 ^ 

n „ 

= n\e n (g) ]T sgn(i/) J] J dxdy(l - x)*-^ 1 ^ - y^-^e-^ 

where we have set z/ = rcr -1 , with the same signature as ar, and explicitly factored out 
the sum over a. Moreover, the dummy integration variables have been rebaptized x and 
y. The partition function takes therefore the form 



Z n (a,b) = nle n (g) <p n (a,b) 
where (f> n (a, b) is the n x n determinant 



(3.4) 



(j) n (a,b) = det 


J dxdy(l - xy-^il - y)3-b-l e -Ngxy 








l<i,j<n 



(3.5) 
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3.2. Discrete Hirota equation 

The technique used in this section to derive a discrete Hirota equation for cf> n is 
borrowed from [12|j . The argument is based on a general quadratic relation satisfied by 
the minors of the determinant D of any (n + 1) x (n + 1) matrix. If we denote by Dij the 
n x n minor of D obtained by erasing the i-th row and j-th column, and Di 1 ^ 2 .j 1 j 2 the 
(n — 1) x (n — 1) minor obtained by removing the rows and columns ji,j2, we have 
the relation 

-D -Dl,n+l;l,ra+l = ^n+1,71+1 -^1,1 — -^l,n+l -Dn+1,1 (3-6) 

This is proved for instance by collecting the monomials in the expansions of the various 
minors. This may also be viewed as a particular case of the Pliicker relations ||13|| . When 



applied to (p n+ i(a + 1,6+1), (|3.6|) immediately translates into the quadratic relation 

n+ i(a+l,6 + l)0 n _i(a,6) = n (o + 1,6+ l)0 n (a,6) - 0»(a,6+ l)0 n (o + 1,6) (3.7) 

where we have absorbed the various restrictions on the matrix indices z, j into correspond- 
ing shifts of a and 6, while the index of refers to the size of the minor. The equation 
(|3.7|) , is known as a discrete Hirota equation, playing a central role in integrable systems 
(see |l3j for the study of analogous equations). 

We may recast (|3.7|) in terms of Z n , using ( |3.5|) and the value ( |2.49| ) of e n (g) satisfying 

(n + l)\e n+1 {g)(n - l)!e w _i(ff) _ 1 n\ _n_ . , 

(n!e n Q7)) 2 Ng{n-l)\ iV^ 1 ' ] 

This leads to 



n— Z n+1 (a + 1,6+1) Z n _i(a,6) 

= Z n (a + 1,6+1) Z n (a y b) - Z n (a,b+ 1) Z n (a + 1,6) 

' (3.9) 

The equation (|3.9|) gives a completely different information on Z n (a,b) or f n (a,b) than 
the recursion relation (|2.52j ). Indeed, the latter determines the dependence of Z n on a, 
6 from the initial data Z\ and through a recursive process, the former instead gives a 
three-dimensional recursion relation on the parameters a, 6 and n. In terms of the free 
energy 

f n (a,b) = Log Z n (a, 6) (3.10) 
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and using the difference operator 

o x f(x) 

we can rewrite ( ^.9| ) as 



= f( x + l)-f( x ) 



(3.11) 



S a S b fn(a,b) = -Log(l - n-eSMa+hb+D-s^^iaS) 



(3.12) 



Substituting for / n (a, b) the small t expansion ( |2.60| ), this gives yet another recursion 
relation for the coefficients w n ,fc(o,&), with the form 



$ a SbVn,k+i = <p(wj,i\j e {n + l,n,n-l}; I < k) 



(3.13) 



This is a recursion over k, with the initial conditions w nj o = and w nj i = nab. However, 
it involves a discrete integration wrt a and b, namely we have to solve an equation of the 
form S a Sbl = h at each step. As h is a polynomial of a and 5, this is readily done by 
noticing that the polynomial u n k(a, b) is always a factor of ab. Indeed, from the original 
interpretation of f n =N as a sum over connected fatgraphs as the logarithm of ( p..5|) , each 
connected fatgraph must have at least one M and one -R-vertex, resulting in a factor ab. 
For n 7^ N, an analogous expansion holds, and the same conclusion applies. In this way, 
we recover the results ( |2.63| ) . Note that the step of discrete integration is simpler here than 
in (|3.13|), where the two discrete integrations wrt n actually involved the determination 



of u n k for all n = 0, 1, 2, 



Here this step can be done in general without having to 



compute explicitly for all (integer) values a, b = 0, 1, 2, ... It is much more economical when 
implemented on a computer. 



3. 3. Large N limit 

In the planar iV —>■ oo limit, we still set z = n/N. We also assume the leading behavior 



f n ~ iV 2 f(z;p, q,t), so that the large iV free energy ( |2.65| ) reads Fo(p,q,t) = f(l;p,q,t). 
The above difference operators ( ^.11| ) become differential operators, namely 



5 n — > —d z 



(3.14) 
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and the Hirota equation Q3.12 ) becomes a partial differential equation for the function 



d p d q f = -Log(l - tze d '( d ' +d * +a J f ) 



(3.15) 



Note that / is not determined entirely by this equation, in particular we may add to 
/ any function of t. However, as mentioned earlier, / has pq in factor, which implies 
that d p f(z; 0, 0, t) = d q f(z;0,0,t) = 0. We therefore solve an equation of the form 
dpd q f(p,q) = g{p,q) as / = L dx f{j dy g(x,y) (we will actually implement this on the 
coefficients Uk(z;p,q) of the expansion (|2.6U[), which are polynomials of p, q for all k). 



Substituting the small t expansion (|2.71|) into ( |3.15|) , and integrating wrt p and g, we 
finally obtain 

r Wn (z; M )f = - f dx rrfyLog(l-^-o az( ^ + ^ +a! ' )wfe(z;x ' y)tfe ) (3.16) 

n>0 J ° J ° 

which is a recursion relation for io n . With the initial condition ujq(z;p, q) = 0, we recover 
the result ( |2.84| ). Note that both ( |2.68|) and ( |3.15| ) are partial differential equations deter- 
mining / up to initial conditions, they must therefore be compatible, a highly non-trivial 
fact considering their very different forms. 

Remarkably, combining both equations (|2.68|) and ( |3.15| ) allows us to solve the case 
p = q = z explicitly, with the following result for the reduced free energj 




as we will prove now. The (p, q, z) symmetry Q2.86| ) of / implies that 



d z (d z + d p + d q )f(z;z,z,t) = (d 2 z +2d p d q )f(z;z,z,t) = ^d 2 z F(z,t) 



(3.17) 



(3.18) 



Now using the equations (|2.68|) and (|3.15| ) in the form 

t 

z 



d 2 j = Lo g (i + -(ta t ) 2 /) 

d p d q f = -Log(l- zte d ^ d * +d *> +d ^ f ) 



(3.19) 



This result coincides with the generating function of planar rooted bicubic maps of In 



vertices [14]. 
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and adding twice the second line to the first one, we get, using ( |3.18| ), the equation satisfied 
by F(z,t) 

\d 2 z F(z,t) = Log(l + t -(td t ) 2 F(z,t)) 

z (6.20) 

-2Log(l -zteWW) 



To check that (|3.17| ) is a solution of this equation, we can relate the expression 

ld 2 z F(z,t) = |>z) fc ^^ (3.21) 
to the generating function of the celebrated Catalan numbers 

C(x) = = J2c k x k (3.22) 

with 

(2*0! r , 9 ,> 

Ck = WTTy. (3 - 23) 

Setting x = 2zt, we may reexpress the result (|3.21| ) as 

o 2y \y/l-4y J (3.24) 

2 \1 + VI - 4x / 2 
= LogC(x) 

Exponentiating both sides of (|3.20|) and using ( p.24|) , our result will hold provided 

C(x)(l-^C(x)) 2 = l + -(tdt) 2 F(z,t) (3.25) 

with x = 2zt as above. But thanks to the quadratic relation xC(x) 2 = C(x) — 1 satisfied 
by (|3.22|) , we may rewrite the lhs of ( |3.25| ) as 

1 + ^C{xf = 1 + jC(x)(C(x)-1) 
= l + i((i-x)C(x)-l) 



x k 



1 + ^2~^( c k -Cfc-l) 
fe=l 

l + A^ tz) (*-!)!(* + 2)1 

i + 3tz ( tao 2 f;(2tz)^^-^ 



(3.26) 



£;!(£; + 2) 
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This proves (|3.25|) for F(z,t) as in ( j3.17|) . Since ( |3.17| ) also satisfies the correct initial 



conditions F(0,t) = /(0;0,0,t) = and d z F(0,t) = 3d z f(0; 0, 0, t) = 0, it is therefore 
identified as the desired solution for p = q = z. Note that, from the convergence radius of 



the series (|3.17| ), this solution has a critical point at 



U{z,z,z) = l- (3.27) 
8z 

Actually, the generating function &tp = q = z = l for bi-colored fatgraphs with a 
marked edge reads (see Appendix A for a similar study) 

oo 

(td t )F(l,t) = J2 tk »k (3.28) 



k=i 



where 



Remarkably, the numbers have already emerged in the context of arches and meanders 
|I~6|1 , within the framework of the Temperley-Lieb algebra. A more detailed study would 



indeed lead to the identification between certain pairs of arch configurations of order k (with 
k arches) and bi-colored diagrams with k edges, thus yielding a different interpretation of 
the result ( ggg) . 



4. Direct expansion 

4-1. The partition function as a sum over integers 

We start from the formula ( |3.5|) for the n x n determinant <f> n (a,b). For generic 
(non-integer) values of a and 6, we may expand 

(1 _ x) i-a-i _ ^ T(k + a-i + l) x k 

1 J ~ 2^ r(l + a-z) k\ [ ' 

k>0 v ' 

and similarly for (1 — . Using the formal integral 

/ ' dxdyx a y e~ Nxy/t = a\ 8 a j (t/N) a+1 (4.2) 
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for any integer a, we get 



(a, b) = det 



1 T(k + a - i + l)T(k + b - j + 1) ( t , k +i 
L ^fc! T(l + a-z) T(l + b-j) 'iV' 



l<i,j'<n 



_ ^ 1 m + a-i + 1) (t/N)^ (43) 

- 2. life., r(i + a-t) r(i + 6 - o den {kl + b 3 + 1) i^,j<n 

where we have used the multilinearity of the determinant to extract line by line the sum- 
mations over fc's. Factoring T(ki + b — n + 1) out of each line (number i) of the remaining 
determinant, we are left with a determinant of the form 

det [(k i + b-j)(k i + b-j-l)...(k i + b-n + l)] = det [q n -j(ki)] (4.4) 



where the polynomials q m {x) = x m +lower degree are monic. As in ( |2.5| ), we may identify 
(fO]) as the Vandermonde determinant A(/c) of the diagonal matrix k =diag(/ci, /c2, k n ). 
This yields 

, / ,\ T T/# /-\T\k-A-i 1 + a - z + 1) r(fci + b - n + 1) , 

A.M) = E n(^) fa+1 ^ V(i + a-o r(i + t-i) AW (45) 

fcl ,...,k n >0 i=l 

Using the antisymmetry of A(fc), we may symmetrize the above expression into 
A, I M - I V" FT 1 (t/iy) tl+1 r(fc i + 6-n + l) 

W ' ' = "■ h ..^ a „iiwr(i + -0 r(i + »-i) (4 . 6) 

x A(k) det[T(ki + a-j + l)]i<ij< n 

Factoring again F(ki + a — n + 1) for each line (number i) of the last determinant, and 
repeating the above trick, we finally get 

mm) = i e ^n(^^ r( ^;:_",; 1)r( ^:;_",; 1) ^ 

fci,...,fc„>0 i=l v ; V 1 ; 

or, using (|3.4j), 
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The expression summed over is symmetric in the fcj's, hence we may further reduce the 
sum to only strictly increasing sequences of fcj's, namely 



Z n (a,b) = 

(t/N)^-^ 1 F(ki + a-n + l)T(ki + b-n + l) 



e A ^ 2 n 



ki\ T(l + a-i) T(l + b 

0<k 1 <k 2 ..-<k n i=l v ' K ' K 

(4.9) 

This expression leads to a straightforward expansion of Z n in powers of t. 

The leading t term corresponds to ki = i — 1, which we refer to as the fundamental 
(or groundstate) configuration of the fcj's. It has the contribution 1. The term of order one 
in t is obtained by creating the excited state where ki are unchanged for % = 1, ...,n — 1 
and k n = n — 1 — > k n = n. Its contribution is easily computed as nabt/N. We may think 
of the power of t as the energy of the excited state (1 here). In general the term of order t k 
will be obtained by a number of excitations of the form ki = i — 1 + Ski, which respect the 
strict ordering of the k's (i.e., Ski < Ski+i for all i), and with a total energy k = J2i 
This makes the expansion of Z n (a, b) quite explicit. 

4-2. Explicit Expansions 

As an illustration of the use of (|4.9[), let us re- work the example a = — 1, b arbitrary, 



of (|2.56|) , by using the expansion (|4.9| ). The terms T(l + a — i) = T(—i), % = 1, 2, n, in 
the denominator may cause the contribution to vanish, unless they are counterbalanced 
by terms of the form T[—j), j > coming from the numerator. Hence, the only non-zero 
contributions to the sum (|4.9|) arise when ki + a — n + 1 = ki — n < 0, for all i. The k's 
being strictly ordered, they are either in the fundamental configuration 

K = (fci,fc2,fc 3 , ...,*») = (0,l,2,3,...,n-l) (4.10) 

or in one of the following n — 1 excited states 

Kj = (0,l,...,j-2,j-l,j + l,j + 2,...,n) (4.11) 

for j = 1, 2, Ti—l. Assembling the contributions from the configurations Kq, Ki, K n -\, 
we get 

26 



which is equivalent to ( 2.55 ) . 

The example a = 1 is also instructive. The denominators T(2 — i) cause trouble for 
i = 2, 3, n. Hence one of the /c's may be chosen freely, while all the others must satisfy 
ki + a — n + 1 = hi — (n — 2) < 0. The strict ordering of the /c's imposes the choice 

(k 1 ,k 2 , ...,k n -x) = (0,1, 2,..., 7i- 2) and k n > n - 1 (4.13) 

This results in an infinity of possible excitations k n = n — 1 + k above the fundamental 
k = 0. These contribute for 

which is equivalent to (|2.58| ) . 

The formula ((OJ) permits also to prove directly the previously-mentioned (a, b, n) 
symmetry of Z n (a,b) = Z a (n,b) for a = m any positive integer. Indeed, for a = 1, if we 
compute Zi(n, b) using (fO|), we find that k\ = k may take the fundamental value k = 
or any excitation k > 0, resulting in the same sum as ( [4.14j ), therefore 

Zi(n,6) = Z n (l,6) (4.15) 



More generally, let us evaluate Z n {m 1 b) using (|4.9|) . The term T(l + a — i) = T(m + 

in the denominator causes problems for i > m + 1 only. Assuming that n > m, this forces 

(n — m) of the /c's to satisfy ki + a — n + 1 = ki — (n — m — 1) < 0. These are necessarily 

(fci,fc2, — ,k n -m) = (0, 1, n - m - 1) (4.16) 

as the sequence is strictly increasing. The remaining /c's may be in either the fundamental 
state or any excited state above it, namely 

< ... < k n (4.17) 

If we look instead at Z m (n, b), there are no bad denominators, and all the excitations of 
the m /c"s are permitted, namely any 

< k[ < k' 2 < ... < k' m (4.18) 

The two sets of excitations ( f4.17j ) and ( [1.18|) can be mapped onto each other by writing 

kj = n — m + kj +m _ n V j = n — m + 1, n — m + 2, n (4.19) 



and it is a straightforward exercise to check that their contributions in ( ^.9| ) are identical. 
This proves the symmetry Z n (m, b) = Z m (n, b) for all m, n, b. 
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4-3. The large N limit: general solution 

When N is large, with n = zN, we expect the sum ( |4.8|) to be dominated by integers 
of the form fcj = Ncti, where ai are real numbers. The index i itself ranging from 1 to 
n = zN, we may let it scale as i = Ns, s ranging from to z. The sum (|4.8[ ) may therefore 
be approximated by the multiple integral 



Z n (a,b) ~ N n da 1 ...da n e N2s{ai '- a ^ q ^ 
Jo 

where we have identified the effective action 

1 n 

S(a>i, a n ;p, q, z) = — ^ [( a i + V ~ ^)(Log(a; + p - z) - 1) 



(4.20) 



i=l 



+ (cti + q- z)(Log(ai + q - z) - 1) - a^Log^) - 1) + a^Logt] 



+ 



72 E Lo ^ 



N 2 



\oi4 — a. 



/ ds[(p- s)(Log(p- s) - 1) 



(4.21) 



+ (g - s) (Log(<? - s) - 1) + s(Log(s) - 1) + s Log t] 

by use of the Stirling formula. The integral ( 4.20|) is dominated by the saddle-point 
dS/dcti = 0, namely 

(ai + p - z)(aLi + q - z) \ 2 1 

/V J—~< rv, — rv„- 



Log t 



a.; 



(4.22) 



for z = l,2,...,n. Introducing the large N density 

1 n 

iV — >oe 'v * — * 

with the normalization 



i=l 



J dap(a) = z 



(4.23) 



(4.24) 



we get the limiting equation 



Log t 



(a + p — z)(a + q — z) 



2P / dp 



m 

13 — a 



(4.25) 



where the symbol P stands for the principal value. As usual, we introduce the resolvent 



d(3 



a — (3 



(4.26) 
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for any complex number a. In particular, u{a) behaves for large \a\ as 

u(a) ~ - / p((3)d(3 = - (4.27) 
a J a 

by use of the normalization (|4.24| ). Moreover, the saddle-point equation Q4.25|) may be 
written as 

/ \ / \ ( (a + p — z)(a + q — z)\ , , 

lim u(a + ie) +w(a-ie) = -Log ^ '- (4.28) 



If we can solve this equation, we can deduce the large N value of the free energy from 



the effective action (|4.21| ). Indeed, writing S(a*;p,q, z;t) = f(p,q,z;t) the saddle-point 
action, we may compute 

td t f = £d«^)^U =a *+^ 

' 2 (4-29) 

= J /3p(J3)dp - Z - 

where we have used ( |4.22|) at the saddle-point ati = a* and the definition ([4.23Q of the 
limiting density p. In terms of the resolvent 

00 O 

u(a) = J2 ~, (4-30) 



cr 

1=1 



this simply reads 



z 2 



td t f(p,q,z;t) = n 2 - — (4.31) 



To solve ( 4.28 ), let us perform the change of variables 

a = z — r — 25cosh </> 

p = r + 25cosh cf>i 

q = r + 25cosh <p 2 

z = r + 25cosh <p 3 



(4.32) 



where r and 6 are two free parameters, to be fixed later. We will take the variable (p = i(p, 
G [0, 7r], to be purely imaginary on the support of the distribution p, so that z — r and 
45 measure the center and width of this support. We define 

u = e~*, T = tanh(0/2) = ^—^ 

1 + U l-u- ^ 
Ui = e-* 4 , T t = tanh(0 i /2) = — —, i = 1,2,3 

1+Ui 
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Then we have 

L ia + p - z)ia + q - z) = 6tj J^0^LzI^0_ (4.34) 

Since on the support of p, T = itan(<£>/2) satisfies T* = — T, the saddle point equation 
( 4.28Q has the following solution 



2(T 1 + T)(T 2 +T)(1-T 3 ) 



for the resolvent provided we take 5 such that 



(i-Ti)(i-r 2 )(i-r 3 ) (a 

ot = -, -—r, ^ttt ^TT = uiu^us (4. do) 

(1+Ti)(l+T 2 )(l+T 3 ) 1 J 

Moreover, the large a behavior ( |4.27| ) imposes that 

I xu 

u(a) ~ w(wi + U2 ) = — — (4.37) 

u 3 d 

namely that 

uiU2U 3 ( u\ — U2) = zt (4.38) 

U3 

Using the third line of ( |4.32| ) , this fixes the value of r to 

r = -8{ux +U2 + U3) (4.39) 



Substituting this value into the second and third lines of (|4.32|) , and using the value (|4.36Q 
for 5, we get 



/ 1 

UiU 2 U 3 ( U2-U3) = pt 

Ul 

/ 1 

uiu 2 u 3 ( ui - u 3 ) = qt 

u 2 



(4.40) 



Introducing the new variables 

Ui = u 2 u 3 U 2 = U1U3 U 3 = U1U2 (4-41) 
we may recast ( |4.38|) and ( |4.40|) into 

i — U 2 — Us) = pt 

(4.42) 



EMI 


-u 2 


-U 3 ) 


= pt 


u 2 (i 


-Ui 


-U3) 


= qt 


u 3 (i 


-Ui 


-u 2 ) 


= zt 
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By eliminating U 2 and U 3 as 



U2 ~ 2 (1 "^ + T^T- 



1(1-^ 
2^ C/i 



(4.43) 



this reduces to the following 5-th order equation for U\ = £/i(p, </, z; t) 



U*(l-U 1 ) 2 (l-2U 1 + 2(p-q-z)t) = t 2 {{l-U 1 ) 2 p 2 -Ul(z-qf) 

(AAA) 

We have to retain the unique solution behaving for small t as U\ ~ pi (this is the unique 
small t behavior compatible with Eq.(4.49) below). The other U's are obtained by ex- 
changing the values of p, q, z, namely 

U 2 (p,q,z;t) = Ui(q,p,z;t) 

(4.45) 

U 3 {p,q,z;t) = Ux(z,q,p;t) 
We then have to plug these solutions back into 



5 = jVUiU 2 U 3 

r = -\{U X U 2 + U 2 U 3 + U l U 3 ) 



(4.46) 



U; 



VU1U2U3 
Ui 



1,2,3 



For instance, the density associated to the resolvent Q4.35Q reads 

P( a ) = TT^ - nm u ( a + ~~ u ( a ~ 

2n\ e^o+ 

= -J-(^-2(*i + * 2 + # 3 )) 
Ztv 



(4.47) 



with a = z — r — 2<5cos(<£>), and we have introduced the quantities 



(4.48) 



The free energy is easily obtained from ( f4.31|) by expanding the resolvent (^4.35| ) up to 
the second order in 1/a. We find 



(4.49) 




Note that this is explicitly symmetric in p, q, z as expected. 
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4-4- Particular cases 



In the passage ( |L43|) from (|4.42|) to ( 4.44T) , we implicitly assumed that U\ was different 



from and 1. Let us examine these two cases in more details. 

If U\ = 0, we read from ( |4.42|) that p = 0. Conversely, if p = 0, we deduce from Q4.42| ) 
that either U\ = identically or U 2 + Us = 1 identically. The latter solution is however 
not compatible with the requirement that U 2 + Us ~ (q + z)t at small t and should be 
discarded. Setting U\ = in ( |4.42| ) leads to 

U 2 - U 2 U 3 = qt 

(4.50) 

Us - U 2 U 3 = zt 

Since 5 = in this limit, the support of the distribution p reduces to one point located at 

a = — - — (4-51) 



which, using ( |4.50| ), is solution of 

t(a + q)(a + z) = a (4.52) 

We thus recover the saddle-point result (|2.75|) with p <-» z. 

If Ui = 1, we read from ( [4.42| ) that q = z. Conversely, if q = z, we deduce from ( ^4.42| ) 



that either U\ = 1 identically or U 2 = Us identically. From the requirement that U± ~ pt 
for small i, we see that it is now the latter solution which has to be kept while the solution 
U\ = 1 is simply unphysical. Let us nevertheless study this particular case q = z, now 
described by its correct solution U 2 = U3. We now get a third order equation for U\ 

Q(U 1: t) = Uf(l-2U 1 + 2(p-2z)t) -p 2 t 2 = (4.53) 

This fixes Ui as a function of t. Writing d t Ui = —d t Q(Ui,t)/d u Q(Ui,t), we find a first 
order critical point t = t*(p,z,z), Ui = m*, whenever d u Q(u*,t*) = 0, hence 

d u Q(u*,U) = 2m*(1 -2«* + 2(p-2z)£*) -2ul = (4.54) 

By combining with ( ^4.53j ), we get ul = (pt*) 2 . As > 0, it is convenient to introduce the 
parametrization 

u * = 1 Te ( 4 - 55 ) 



3 
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where 

9 > when - < u* 
4 ~~ 

9 = icp, if > when < u* < — 



(4.56) 



The critical point is then 



t*(p,z,z) = -j (4.57) 

Hp cos | 



with e = ±1. The values of 9 and e are fixed by the ratio 

z q 1 + e cos 9 
p p 2 

cos 2 | if e = 1 
sin 2 f if e = -1 



(4.58) 



For z/p > 1, ( |4.58| ) has exactly one solution corresponding to e = 1, z/p = cosh 2 (<p/2) 
and we find a critical point at £* = l/(8pcosh 3 (ip/3)). We thus have 

U(p, z, z) = -j— (4.59) 

8p cosh ( | arccosh . ~ 



For z/p < 0, C |4-58|) gives also exactly one solution (u*,t*) with now e = —1. However, this 



point is not on the correct branch of the solution of ( f4.53| ), i.e. that satisfying U\ ~ pt at 
small t. This is easily checked by solving explicitly ( f4.53| ), which is quadratic in t in the 
form t(Ui) rather than Ui(t), and choosing the solution with the correct small t behavior. 
We thus deduce that there is no critical point (R = +oo) for z/p < 0. 

For < z/p < 1 (|4.53|) has an infinite number of solutions of the form e = 1, 9 = 



# + 2/c7r, 9 = (2k + l)Tv-9 and e = -1, 9 = # + (2/c + l)7r, = (2k + 2)7i-9 with ^ e [0, tt[ 
and fc a non negative integer. This leads to exactly three candidates for the critical point 
t« = l/(8pcos 3 (^ /3)), t* = l/(8pcos 3 (^ /3 + 2tt/3)) and = l/(8pcos 3 (^ /3 + 4tt/3)). 
Here again, it is easy to check that the critical point attached to the correct branch is the 
first one, namely 

U (p, z, z) = j- (4.60) 

8p cos 3 f g arccos J | J 

where the arccos is taken in the range [0, 7r[. 

When p = q = z, we recover the critical value £*(z, 2) = l/8z of (|3.27|) by setting 



9q = 0, with moreover it* = 1/4. This case p = q = z can be solved more directly by 
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noting that, from ( |4. 42| ) , we then have U\ = XJ% = Us = U, where U is the solution of a 
quadratic equation £7(1 — 2£7) = zt, namely 

£7 = 1(1-^1-82*) = ztC(2zt) (4.61) 

The corresponding free energy reads 



£7 3 

td t F(z;t) = — (l-3£7) = — ((8zt - l)C(2zt) + 1 - 6zt) (4.62) 



which coincides with our previous result ( |3. 28| ) ( pT29| ) . 



4-5. Critical Points 



The solution U\ to (|4.44|) such that U\ ~ pt when t — > will in general have a 
convergent series expansion until we reach a critical point t = t* (p, , corresponding to 
the convergence radius R = t*. To determine this point, we write the (finite) Taylor series 
of the polynomial 

P(u, t) = u 2 (l- u) 2 (l - 2u + 2(p - q - z)t) - (1 - u) 2 p 2 t 2 + u 2 (q- z) 2 t 2 (4.63) 
namely 

^ \^ 

P(u,t) = P{u*,U) + (u - u*)d u P{u*,U) + K - -^-d%P{u*,U) . . 

2, (4.64) 

+ (t-U)dtP(u*,U) + ... 
Imposing that P vanishes, we get a first order critical point 

~ (t*-t) 1/2 (4.65) 

by demanding that 

P(u*,U) = d u P{u*,U) = (4.66) 

while 

d 2 u P(u„ U) ^ d t P(u„ U) ^ (4.67) 
This means that the polynomial P(u,t*) takes the form 

P(u,U) = (u-u*) 2 (-2u 3 + au 2 + bu + c) (4.68) 
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for some numbers a, 6, c to be determined. Introducing the parameters 
a = 2(p - q - z)U P = P 2 tl 7 = (q ~ z)Hl 
and identifying (|4.63|) with (|4.68|) , we get 

a + Au* = a + 5 
b - 2au* - 2ul = -2(a + 2) 
c — 2bu* + aul = 1 + a — /3 + 7 



(4.69) 



(4.70) 



bu* — 2cu„ 



cu* 



2/3 
-P 



Eliminating P from the last two lines of ( |4.70| ), and a from the first two, we find 

1 — w* 



2c- 
3 



c 



(4.71) 



Here we can assume that it* 7^ 0, 1, since the situation where U\ = or 1 has been treated 
in the previous section. The equations ( f4.71|) translate into the following parametrization 
of a, P, 7 of (|4T69|) 

a = 2(p — q — z)t* = 3it* — 2 + £ 

= P 2 tl = £ul (4.72) 

7 = (q~z) 2 tl = (l-ti^Cl-O 
where we have set c = — This may be viewed as a parametric curve (pt*, qt*, zt*) = 
h(u*, £), where and £ are real parameters, with > and (1 — — £) > 0. This 

in turn permits to express for instance 



(4.73) 




where e =sgn(p£*), as an implicit function of the ratios q/p and z/p, where 



q 
p 

z 
P 



2 + ^ 



1/^ ^,1 — /(l — — e 3u 5 
2^ «■* V w*£ 2u*y/£u 



2 + ^ 



1/ / (1 - - 3u 5 

2^ w* V w *£ 2u^y/^u 



(4.74) 
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and where e' =sgn(p(<j» — z)). Note that the transformation (£, u*) — > (1 — £, 1 — m*) simply 
amounts to 

p -> ^ e'(g-z) 

9 -> ^p+^Q-e'^ (4.75) 

z — ► + — e '-2 



We must now check the validity of this critical point, by making sure that (|4.67|) is 
satisfied. At the critical point, we find 

d t P(u m ,U) = j-ul(l - «*) 2 (it* - £) 

** (4.76) 

d u P{u*,t*) = 6u*(l — it*) (it* — £) 



Provided we take it* 7^ 0, we see that ( |4.67|) will hold. As mentioned above, the case 



u* = corresponds to p = and has already been studied in the previous section. Taking 
it* = 1 implies that g = z, a particular case also treated in the previous section. Note 
that in this case, the correct solution for the critical point is not at = 1, which lies 
on a wrong branch, but at some < 1 given by ( |4.55| ). The case 11* = £ corresponds 



to symmetric counterparts of these particular cases. Indeed, from (|4.72Q , we find that 
pt* = eu 2 , (q — z)t* = ee'(l — u*) 2 and (q + z)t* = eu 2 + 1 — 2ii*, hence 

(qU,zQ = ((l-w,) 2 ,0), (0,(l-w*) 2 ), (1-2^,-n 2 ), (-n 2 ,l-2^) (4.77) 

according to the values of the signs (e, e') = (+, +), (+, — ), (— , +), (— , — ). The first two 
cases correspond to z = or q = and are equivalent to the p = particular case under 
the change p <-> z or p <-> The two last cases correspond to either p = q or p = z, which 
are now equivalent to the q = z solution. Here again, it* = £ lies on a wrong branch of the 
solution. 

From now on, we shall ignore the particular cases by demanding that none of the 
parameters p, q and z vanishes and that they are all different. Using the (p,q,z) symmetry, 
we can furthermore assume without loss of generality that —1 < z/p < q/p < 1, with in 
particular e' = 1. For fixed q/p and z/p in this range, any solution of Q4.74|) is a candidate 
to describe a first order critical point, of the form ( |4.65| ), with t* given by ( |4.73| ). We still 



have to make sure that both terms in (|4.67|) have the same sign, in order for the solution 



to read 11* — u ~ (1 — t/t*) 1 / 2 . This is the case if and only if 

< < 1 (4.78) 
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which in turn implies < £ < 1. The parametrization Q4.74j ) is further restricted by 
demanding that the critical value u* indeed corresponds to a solution U\(t) = pt + 0(t 2 ) 
for small t. Using the equation (|4.44j) , we may solve for t as a function of U\. We have to 
pick the right branch of this quadratic equation, namely that for which t{U\) = ^y+0(U 2 ) 
for small U\. It reads 



pt(u) 



u(l — u) 



(1 - „)2 _ (az«) u< 



1 u(l — u) 

p 



+ 



(4.79) 



u 2 (l - u) 2 1 - 



q + z 
p 



+ (l-2u) 



q - z , 2 

p 



In this formulation, a critical point corresponds to a local extremum of pt(u). The 
parametrization ( [4.74 ) corresponds to this branch provided 



sgn(£ - 



(4.80) 



This is obtained by substituting (|4.74j) into ( |4.79|) , and comparing the result with ( |4.73|) . 
So finally, we have fixed the values of e and e' in ( |4.74[ ). To ensure the uniqueness of the 
critical point we found, we have to check that the parametrization ( [4.74j ) is bijective. We 
must therefore compute the Jacobian J of the map (it*, £) — > (q/z,p/z) induced by (f4.74|) : 
we find 

Z . /i _ u 

(4.81) 



J 



This is non- vanishing for < « ^ < 1. We therefore have for fixed q/p and z/p at 
most two acceptable solutions of ( |4.74j ), one in the range < < £ < 1 with e = 1 and 
one in the range < £ < u* < 1 with e = — 1. This latter solution can be eliminated 
by noting that it corresponds to a value of £* given by ( f4.73|) such that pt* < 0, which is 
unphysical. Indeed, starting from pt ~ U\ at small t, one cannot reach continuously an 
extremum (u*,t*) of t(u) with u* > and pt* < without first reaching an extremum 
at some it* > 0, pt* > 0. This fixes e = 1 and < u* < £ < 1, which leaves us to 
at most one solution of (|4.74|) . In this sense, the parametrization ( [4.73|)(|4.74|) leads to 
a unique expression t#(p,q,z) = (l/p)T(q/p, z/p) for those values of q/z and z/p in the 
range — 1 < z/p < q/p < 1 for which a critical point actually exists. 

To determine which values of q/p and z/p give rise to a critical point, we distinguish 
two situations in the formula ( |4.79| ), as described in more details in Appendix B. For 
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z/p < 0, the discriminant in ( 4.79 ) does not vanish. The function pt(u) is a strictly 
increasing function from i = at w = to t = +oo at the value u = 1/(1 + (q — z)/p) 
where the denominator vanishes. There is no critical point in this case, hence R = +oo. 
For z/p > 0, the discriminant in ( |4.79| ) vanishes at some value uq of u strictly less than 
the limiting value 1/(1 + (q — z) /p) for the vanishing of the denominator. At uq, the slope 
of t{u) is — oo, which implies the existence of a critical point at some between and uq 
with pt* > 0. 

In conclusion, we find a critical point only if < z/p < q/p < 1, i.e. if p, q and z have 
the same sign, which in turn is also the sign of £*. 



z 
P 



8»cos 3 ( — arccos 

V3 Vp 




S-zcosh' 3 I — arccoshW — 



1 
P 



Fig. 2: Critical value t*(p,q,z) in the domain —1 < z/p < q/p < 1. A 
critical point is reached only in the sub-domain z/p > (shaded region). The 
limiting values of t* are explicited at the boundary of this sub-domain. 

Our results are summarized in Fig. 2, where we have extended the domain to — 1 < 
z/p < q/p < 1 by reintroducing the results of the previous section for the particular cases. 

Using the results ( fl.76| ), we may finally compute the first sub- leading coefficient in 
the expansion 

Ui = u* -k(U (4.82) 
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namely 



K 



'it*(l — u„ 

3^ 5k 



(4.83) 



The free energy follows from ( f4.49|) . Actually, it is simpler to compute 

t 2 (td t ) 2 f = U X U 2 U 3 (4.84) 
easily obtained by use of Q142|) . Substituting the values Q4.43|) for {7 2 and C/3, this reads 



*w/ = £((i-£) a M ~ z)t 



l-Ut 



) 2 ) 



.((l + fp-g-^^-pi-^ 



(4.85) 



Substituting the expansion (|4.82|) , we get the expansion 




with a non-vanishing singular term if u* 7^ 0, In particular, for p = q 
u* = 1/4 and f = 1 so that (|4.86|) yields 



(4.86) 
z, we have 



t 2 (td t ) 2 f(z,z,z,t) = — -— VT=m+... (4.87) 

in perfect agreement with (|3.25|) ( |3~2^ ). 

The above result ( [4.86| ) leads to the string susceptibility exponent 7 s t r = —1/2, gov- 
erning the critical behavior of the free energy / ~ (i # — t) 2 ~ lstr . This coincides with the 
string susceptibility exponent of pure gravity, attached to the plain enumeration of planar 
diagrams. This universal result holds for all the values of parameters (p, q, z) where a criti- 
cal point is reached, except in the particular cases where one of these parameters is exactly 
zero, in which case we find the behavior fl2.82p . The convergence radius R = £* may be 
related to the entropy s of three-coloring of the vertices of random triangulations, namely 
R = e s . Actually, as a consequence of (|4.86| ), the contribution to the free energy with a 
fixed number of black- white edges n (= 1/2 x the number of triangular faces) behaves for 
large n as 

e ns(p,q,z) 

u) n (p,q,z) j (4.88) 

«2 



where the three-coloring entropy per black-white edge reads 

s{p,q,z) = -Logt^p^q.z) 

with t* as in (|4.73| ). 



(4.89) 
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4-6. Double Scaling Limit 

We may now reconsider the all-genus expansion of the free energy ( |2.60| ) 

f(p,q,z;t;N) = f n (p,q,l/t,N) = ^ N 2 ~ 2 h <p h (p, q,z; t) 

h>0 



(4.90) 



where the genus h free energy ip^ is a quasi-homogeneous function, satisfying cph(Xp, Xq, Xz; t/X) 
X 2 ~ 2h Lph{p,q, z]t). When t approaches the critical point t^(p 1 q 1 z) of (|4 . 73| ) , the leading 
singular part of (fo behaves as (t* — i) 5//2 , and has a prefactor iV 2 . This suggests to perform 
the double-scaling limit t — > £*, N — > oo, by keeping the parameter 



y 



ATV5 



(4.91) 



fixed. Let us use the expansion parameter e = 1/iV. At the vicinity of t*, let us write the 
total free energy 



f(p,q,z;t;N) = e ' A g (p, q, z) + e &/b g 1 {p,q,z)y + e ' li/>> g2(p,q,z)y 

+ e °/(y; p, 9, *) + e 2/5 ^ 3 (p, <z, ^)y 3 + 0(e 4 / 5 ) 



-2/5 



(4.92) 



This defines the scaling function f(y;p, q, z), capturing the contributions to the free energy 



from all genera in the double scaling limit ( [4.91|) . Its genus zero (y — ► oo) contribution 
reads 

tf \ u *~£ /w*(l - W„) 5 

f{y;p,Q,z) = — - — y y 2 +- (4.93) 

as a consequence of ( |4.86| ). To derive an equation for f(y;p, <?, z), let us start from the re- 
cursion relation (|2.52| ), and expand it in powers of e. With /„ = Log Z n = f(p, q, z; t; N) = 
f(z, t), we find 



f{z + e,t) + f{z-e,t)-2f(z,t) = Log(l + e~ 2 -(td t ) 2 f(z, t)) 



(4.94) 



Substituting the expansion ([4.92[) for /, we finally get at the order e 4 / 5 the equation 



1 

6 L y 



ffif+ifif? = k 2 y 



(4.95) 
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where k is a function of p, q, z only, fixed by the genus zero behavior ( 4.93 ) to be 



k = -(u* -0^3^(1 - w*) (4.96) 

Eq. ( |4.95|) is nothing but the Painleve I equation of two-dimensional pure quantum gravity, 
for the string susceptibility dyf. The dependence on the parameters p, q, z is entirely 
contained in the function k. Note that k is symmetric in p,q,z, as 

k 2 = 75ui 1 \? ) u?\l-vi 1) )(l-u!? ) ){l-u? ) ) (4.97) 

where u* = OS z), = u*(q,p, z), = u*{z, q,p) denote respectively the critical 
values of the functions Ui, U2, Us of (|4.42|) . (Recall that we have written u* as an implicit 
function of p, q, z through the parametrizations (|4.74j ), which must be inverted.) 

We conclude that although our model involves more than the mere enumeration of 
random surfaces, its double-scaling limit satisfies the same universal Painleve I equation 
as pure gravity, up to a simple rescaling involving the function k(p, q, z). 



5. Conclusion 



In this paper, we have investigated the partition function ( 2.64 ) of the vertex tri- 
color ing problem for random triangulations of arbitrary genus. In addition to the exact 
evaluation of the tri-coloring entropy (|4.89|)([4.73 ), we have shown that the model lies in 



the universality class of pure two-dimensional quantum gravity, with a string susceptibility 
exponent ^ str = —1/2. 

5.1. Folding of random triangulations 

It is interesting to note that, in the vertex tri-coloring problem, the coloring itself 
does not give rise to any degree of freedom since for a given triangulation, there is only 
one way to color its vertices (up to global permutations of the three colors) . Varying p, q 
and z only amounts to favoring one of these global permutations. This situation is to be 
contrasted with the tri-coloring of the edges of triangulations. There, even for the regular 
triangular lattice, the number of tri-colorings increases exponentially with the system size. 

Therefore, in our problem of vertex tri-coloring, all the entropy comes from the sum 
over triangulations itself, as for pure gravity. The novelty here comes from the fact that 
this sum is restricted to tri- colorable triangulations. In particular, this tri-colorability 
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constraint implies that all vertices have an even number of surrounding triangles. This 
requirement is clear since the vertices around a vertex of a given color are bi-colored with 
alternating colors. Conversely, this parity condition turns out to be a sufficient condition to 



ensure tri-colorability when applied to triangulations with genus zero, as discussed in JT?]] . 
This is not the case however for triangulations with higher genus. For p = q = z where we 
do not distinguish between the different colors, we have found that the sum over connected 
planar tri-colorable triangulations obeys a quadratic equation and can be given in terms 
of the generating function for Catalan numbers. This result is somewhat surprising since 
Catalan numbers are usually associated with the enumeration of objects which can be 
simply decomposed into two similar smaller objects (like trees, arches, parentheses). We 
do not see any such decomposition here. 

In is interesting to reinterpret our results in the language of folding of the triangu- 
lations. If we imagine a triangulation made of equilateral triangles with unit side length, 
we may want to map its vertices onto those of the regular triangular lattice in the plane, 
with positions ae\ +6e*2 +ce*3, (a, b and c integers) where {e^} is a set of three unit vectors 
in the plane satisfying e i + e*2 + e*3 = 0. Such a mapping will be called a folding of the 
triangulation if it sends nearest neighbors of the triangulation onto nearest neighbors on 
the lattice, thus preserving the equilateral nature of the triangles. Two neighboring trian- 
gles will be either side by side in the plane or on top of each other with a fold between 
them. All triangulations cannot be folded onto the regular planar triangular lattice. The 
possibility for a triangulation to be folded is precisely equivalent to the tri-colorability of 
its vertices. To see this equivalence, we note that the regular triangular lattice can be 
decomposed into three regular triangular sublattices such that any neighboring vertices 
belong to two different sublattices. If we color each of these sublattices with three differ- 
ent colors, this coloring induces a tri-coloring of the vertices of any folded triangulation, 
which is thus tri-colorable. Conversely, if the triangulation is tri-colored, it can clearly 
be mapped onto a single triangle, by sending all vertices of the same color onto the same 
vertex of the triangle. We see here that the tri-coloring of the vertices can be viewed as 
a particular folded configuration of the triangulation where all triangles are mapped onto 
the same triangle. We call this configuration the complete folding since all edges of the 
lattices sustain a fold. Our result simply counts completely folded random triangulations 
with a different weight to the three different image vertices. Clearly, the complete folding 
is only a particular folded state and there are in general many different ways to fold a given 
tri-colorable triangulation, with an image covering in general several triangles in the plane. 
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This corresponds to an additional degree of freedom which is nothing but the tri-coloring 
of the edges of the triangulation. Indeed, any edge of the triangulation is sent onto one 
of the three vectors e±, e*2 or e*3 and we can assign a color to this edge accordingly. In 
conclusion, the problem of folding of random triangulations can be identified with that of 
tri-coloring of both its vertices and its edges, a problem yet to be solved. 

5.2. Hirota equation for multi-matrix models 

The general solutions we found take the form of either some recursion relations ( |2.52| ) , 
( p.9|) or an explicit multiple sum over strictly increasing sequences of integers ( |4.9|) . The 
latter is a discrete analogue of the initial matrix integral itself, as it involves the square of 
the Vandermonde determinant of the integers summed over. Actually the large N limit of 
this sum is dominated by the integral Q4.20| ), which resembles the large iV approximation 
to the n x n one-matrix integral 

" dMe -NU(M;p, q>Z ) (51) 



where the potential U is derived from the effective action S of ( 4.21|) , namely U(M; p, q, z 



NS(Xi, A n ;p, q, z), where Aj denote the eigenvalues of M. This resemblance is however 
not quite an identity, as the corresponding eigenvalue integral should extend over the whole 
real line for ( |5.1| ), whereas it extends only over [0, oo[ in ( |4.20| ). 

The existence of a discrete Hirota bilinear equation for Z n is actually generic for the 
following multi-matrix integral 



Z = J dM 1 dM 2 ...dM k det(l - Mi)-° det(l - M k )~ h e -N™{M u ...,M h ) 



(5.2) 



over n x n Hermitian matrices, with the potential 

k fc-1 



V(M 1 ,...,M k ) = J2 V i( M i) + J2 u i M i M i+ 1 ( 5 - 3 ) 

i=i j=i 



with arbitrary potentials Vi (including possible logarithmic pieces). Indeed, using the 
Itzykson-Zuber integral (|2J]) , we may reduce Z to an integral over the eigenvalues of 
the M^s. Proceeding exactly as in Sect. 4.1, we may recast ( |5.2|) into the determinant 



det 



(5.4) 
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up to an unimportant multiplicative constant. This is clearly a solution of the Hirota 
equation ( |3.7| ). This solution is singled out by its initial value <pi, which is a function of 
the details of the potential V (and in particular of the parameters u±, U2,---Uk-i)- The 
case k = 1 of a one-matrix model is also instructive. In that case, we find that [U2fl 



det 



J dm(l- m y + i- a - 2 e- NV ^ 



satisfies the following bilinear equation 



(j) n +i(a + 2)4> n - 1 {a) = 4> n (a + 2)</> n (a) 



(5.5) 



(5.6) 



In this paper, we have seen however the limits of the use of the Hirota equation, insofar 
as the large iV limit, for instance, is concerned. Because the Hirota equation is only a multi- 
dimensional recursion relation, strongly relying on its initial conditions, it makes it quite 
difficult to extract any asymptotic result. We have had to resort to a different approach to 
reach that goal. Perhaps a suitable mixing of the orthogonal polynomial solution and of 
the Hirota equation could lead to some more definite results: this is still an open question. 

Finally, let us mention that our solution has displayed striking connections to the 
arches and meanders enumeration problems and their reexpressions within the framework 
of the Temperley-Lieb algebra [1^] [16] . We intend to return to these aspects in a later 
publication. 
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Appendix A. The planar free energy for small z 

In this appendix, we compute the leading coefficient (p± (p, q, i) of the free energy 
f(z;p,q,t) when z is small. With a* the solution ( |2.76|) of the saddle-point equation 
Q2.75Q , we have 

¥>i(p,9,t) = S(a*) = pLog(l + — )+ ? Log(l + —)-«* (A.l) 

V Q 



Let us compute ip(p, q, t) = tdt<fi(p, q, t) in terms of a* . Using the formula ( |2.74j) for S(a 



and the fact that, at the saddle point a*, we simply have to take the explicit derivative of 
S(a) wrt t, we obtain immediately 

^(p,q,t) = a* (A.2) 

with a* given by ( |2.76| ). Hence, introducing the reduced variables u = pt and v = qt, we 
have ip(p, q, t) = I(u, v)/t, where the function I(u, v) reads 

I(u,v) = - (l - u - v - y/l - 2(u + v) + (u - v) 2 ) (A.3) 

and satisfies the quadratic equation 

(I + u)(I + v) - I = (A.4) 

Expanding ( |A.3|) in powers of u and v, we find 

k 

I(u,v) = 52Y,IkjV>v k+1 -1 (A.5) 

k>l j=l 

with the integer coefficients 

^if fc V.M (A.6) 



Integrating ifj(p, q, t)/t = I(pt, qt)/t 2 once wrt t, we finally get 



<pi(p,q>t) = ^2t h ti k (p,q) (A.7) 



k>\ 

where 



v*m = lE^y'g^ (Ai 



3=1 
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Fig. 3: Equivalence between (a) a bi-colored fat graph with a unique vertex 
and a marked edge, (b) a system of bi-colored arches and (c) a system of 
arches closed into a set of connected circuits 



and the result ( |2.80|) follows. This result can be checked directly on (|2.84j ). 

From the behavior a* ~ (t* — t) 1 / 2 read on ( |2.76| ), we deduce by integrating (|A.2| ) 
that (pt(p, q,t) ~ (U -t) 3/2 . 

The numbers Ikj ( |A.5|) have already emerged in [15] in the context of arch configura- 



tions. The connection between ifj(p, q, t) and arch configurations can be stated as follows. 
Since (pi (p, q, t) keeps only the linear term in z of f(z, p, q, t) , it is a generating function for 
planar (dual) black and white fatgraphs with exactly 1 vertex. The operation tdt simply 
corresponds to specifying an edge among the E(T) edges of a graph T. We can now stretch 
the unique vertex of the graph into a straight line and let the edges of the graph form a 
system of E(F) arches drawn on top of this line (see Fig. 3). The stretching operation is 
completely specified if we impose that the marked edge becomes the leftmost arch, with 
its white adjacent face above it. The function tfj(p,q,t) is thus a generating function for 
systems of bi-colored arches. Note that the same arch configuration is obtained exactly 
|Aut(r)| times from the different markings of a given graph T, which in turn cancels the 
symmetry factor in the diagrammatic expansion. We thus have 

k 

4(p,q,t) = J2Y, Mk ^ tk P J( i k+1 ~ J ( A - 9 ) 

fc>i j=i 

where Mkj is the number of distinct configurations of k arches drawn on top of an oriented 
line, separating k + 1 domains in exactly j black domains and £ = k + 1 — j white domains 
when colored in black and white with alternating colors, the outermost domain being 
white. 



The number M^j has been computed in |[15|| , with the result 

Mk, 3 = h,j (A. 10) 
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given in ( |A.5|) , with a slightly different interpretation as the number of systems of k arches 
drawn on top of a line which give rise to exactly j connected circuits when closed below the 
line by elementary arches connecting successive points on the line in pairs. This number is 
clearly the same as the number of bi-colored systems with j black domains and k + 1 — j 
white domains by coloring in black the interior of the j connected circuits (see Fig. 3). 

The interpretation of i{j(p : q,t) as the generating function for bi-colored systems of 
arches allows for a more direct computation by use of a simple recursive relation for the 
Mkj. Denoting by 



(A.ll) 



we have the recursion relation 



3=0 



(A.12) 



with the convention Mo(p, q) = q. This equation simply expresses the fact that the leftmost 
arch separates the system of bi-colored arches into two systems of bi-colored arches. Note 
that p and q are exchanged in one of these two systems of arches. The above recursion 
relation translates into a quadratic equation for ^{p, q, t) = J2k>i Mk(p, q)t k 



fp{P,q,t) = t(ip(q,p,t) +p)(ip(p,q,t) + q) 



(A.13) 



which, since ip{p, q, t) is by definition symmetric in the exchange p <-> q 1 is nothing but the 
equation ( |A.4|) . 

Note finally that when p = q, I{u,u) = u(C(u) — 1), where C(u) is the generating 
function of the Catalan numbers ( |3.22| ) ( |3.23| ) , and it follows that 

i! 



(A.14) 



Appendix B. Existence of a critical point 

We start from the solution t±(u) of the equation P(w, t(u)) = 



pt±(u) 



u(l — u) 



(i-ur-( 



1 u(l — u) 

p 



± 



\ 



u 2 (l -u) 2 1 - 



q + z 
p 



+ (l-2u) (1-u) 



q-z , 

p 



(B.l) 
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We want to know whether or not the branch t+(u) (satisfying pt+(u) ~ it at small u) 
reaches continuously by increasing u from a maximum at some value u* between and 
1. We assume here that —1 < z/p < q/p < 1. 

We first note that the denominator in flB.lj ) vanishes for u = u\ with 

u x = (B.2) 

v 



between 1/3 and 1. Let us study the possible vanishing of the discriminant in ( [B.l|) . This 
discriminant vanishes for 

-^MoV^H 2 ^) 2 " 2 ) <■»> 

In the domain < u < 1, the r.h.s. is positive for u between 1/2 and u±, and is maximum 
at 

«2 = 1 - T = (B.4) 

1 + \/¥ 

with the value 

(B.5) 

P J 

The discriminant will thus vanish if and only if 

This is true in the domain —1 < z/p < q/p < lif and only if z/p > 0. In this case, 
the discriminant vanishes exactly once in the range < u < U2 for some value uq with 
moreover 1/2 < u$ < U2 < u\ (since (q — z)/p < 1). At u = Uo, the two branches t±{u) 
meet with an infinite slope (if z/p > 0) and make a loop. Clearly, the branch t+iu) must 
have a maximum at some u* < uq and we find a critical point. For z/p = 0, we find a 
critical point at = u$ = u^- 

Conversely, if z/p < 0, the discriminant does not vanish in the range < u < 1 and, 
since (1 — (q + z)/p) > 0, we find that t + (u) is a strictly increasing function which diverges 
at u = U\. We find no maximum is this case, hence no critical point. 

We give an example of these behaviors in the Figure 4. 
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Fig. 4: The behavior of pt±(u) for < u < 1, here represented for q/p = 2/3 
and decreasing values of z/p. A critical point corresponds to a maximum of 
pt{u) lying on the correct branch of solution, i.e. that leaving the origin with 
a positive slope. For z/p > (fig. a,b,c) this branch forms a loop with the 
wrong branch leaving the origin with a negative slope: a maximum is reached 
on the correct branch before the two branches meet. For z/p = (fig. d), a 
maximum is also reached, lying now exactly at the meeting point of the two 
branches. For z/p < (fig. e,f), the two branches avoid each other and no 
maximum is found for the correct branch. 
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